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ABSTRACT 

We perform time-resolved spectroscopy on the prompt emission in gamma- 
ray bursts (GRBs) and identify a thermal, photospheric component peaking at a 
temperature of a few hundreds keV. This peak does not necessarily coincide with 
the broad band (keV-GeV) power peak. We show that this thermal component 
exhibits a characteristic temporal behavior. We study a sample of 56 long bursts, 
all strong enough to allow time-resolved spectroscopy. We analyze the evolution 
of both the temperature and flux of the thermal component in 49 individual 
time-resolved pulses, for which the temporal coverage is sufficient, and find that 
the temperature is nearly constant during the first few seconds, after which it 
decays as a power law with a sample-averaged index of —0.68. The thermal fiux 
first rises with an averaged power-law index of 0.63 after which it decays with an 
averaged index of —2. The break times are the same to within errors. We find 
that the ratio of the observed to the emergent thermal fiux typically exhibits a 
monotoneous power-law increase during the entire pulse as well as during complex 
bursts. Thermal photons carry a significant fraction (~ 30% to more than 50%) of 
the prompt emission energy (in the observed 25-1900 keV energy band), thereby 
significantly contributing to the high radiative efficiency. Finally, we show here 
that the thermal emission can be used to study the properties of the photosphere, 
hence the physical parameters of the GRB fireball. 

Subject headings: gamma rays: bursts - gamma rays: observations - gamma 
rays: theory - radiation mechanism: thermal 
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Introduction 



Although being studied for several decades now, observations of the prompt emission 
phase of gamma-ray bursts (GRBs) have not yet been able to lead us to a clear understand- 
ing of either the underlying mechanism of the emission, nor its origin. The prompt emission 



as the Band function; 


Band et a 




1993; 


Preece et al. 


1998a. 


2000: 


Barraud et al. 


2003; 


Kaneko et al. 


2006; 


Kaneko et al. 


2008; 


Abdo et al. 


2009; 


Gonzalez et al. 


2009 
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ing phenomenological in nature, in many cases this interpretation was found to be consistent 
with the prediction s of optically-thin synchrotron em i ssion from a power-law distribution of 



energetic electrons (ITayani I Il996l : ICohen et al. 1 119971 : ISchaefer et al. I Il998l : iFrontera et al. 



2000; Wang et al. 20091 ). A peak in the energy spectrum then naturally arises, connected 
to the low-energy cut in the distribution of electrons. Indeed, it is widely accepted that 
the non-thermal spectra arise from the prompt dissipation of a substantial fraction of the 
bulk kinetic ener gy of a relativist ic outflow, originating from a central compact object (for 



review, see, e.g.. 



Meszaros I l2006l ) . The dissipated energy is converted into acceleration of 
electrons, which produce high-energy photons by synchrotron emission and inverse Compton 
(IC) scattering. 

In spite of its many successes, evidence have accumulated in recent years for low-energy 
spectral slopes that are steeper than allowed b y the optically-thin synchrotron or syn- 



chrotron self-Compton (SSC) model predictions (ICrider et al. 1 119971 : iPreece et al. I Il998b 



2OO2I : IChirlanda et al. II2OO3I : IConzalez et"aDl2009h . Moreover, the fittin g is often made to the 



tim e-integrated spectrum. A nalysis of time-resolved spectra done by ICrider et al. I (Il998l ) 
and iGhirlanda et al. I (120031 ) . showed that neither the synchrotron nor the synchrotron self- 
Compton models can explain the instantaneous ORB spectra and their evolution. Motivated 
by these findings, other origins of the emission has been suggested such as small-pitch-angle 
synchrotron emission (Epstein & Petrosian 1973, Med yedev 2006) and inverse Compton scat- 
tering of a self-aborbed seed emission. However, as iBaring fc Braby I (120041 ) have pointed 
out these emission mechanisms rely on almost purely non-thermal distri butions of electrons 
to be able to fit the observed spectra (see also IStern fc PoutanenI ( 120041 )). This is not a re- 
alistic assumption since diffusive shock acceleration typically leads to a strong c ontribution 
of a thermal population. This was re cently shown in the works by, for instance, ISpitkovsky 
( I2OO8I ): ICiannios fc SpitkovskyI (120091 ). In addition, internal shocks are assumed to dissipate 
the bulk kinetic energy, in order for it to be emitted by synchrotro n emission. This proces s 
is however too ineffi c ient to explain th e obse rved energy release (IKobayashi et al.l (119971 ): 
Lazzati et al.l (119991 ): iNysewander et al.l (120081 )) and moreover the predicted correlation be- 
tween peak energy and luminosity, that such a scenario gives rise to, contradicts what is 
observed (Ramirez- Ruiz & Lloyd- Ronning 2002). 
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The difficulties that the purely non-thermal emission models face are readily alle- 
viated by introducing an additional optic ally-thick thermal (blackbody) component that 



may contribute to the qbseryed spectrum (jMeszaros. fc Rees I l2000al : iMeszaros et al. 1 12002 



Daigne &: Mochkovitch 1 120021 : iRees &: Meszaros 1 12005| ) . Two different emission regions con- 
tribute in forming the spectrum as well as its temporal evolution; a narrow thermal com- 
ponent and a broad, non-thermal component are thus combined. In this interpretation, the 
prominent spectral break (peak) in the sub-MeV range is identified with the peak of the 
thermal emission, instead of being related to the radiating electrons at the low-energy cut 
in their energy distributior0. This has significant advantages since it can accommodate (i) 
the existence of hard subpeak spectral slopes, (ii) temporal variations of the spectral shape, 
(iii) the observed spectral correlation between peak energy and luminosity within a burst 
(Golenetski et al. 1983), (iv) the high radiative efficiency required, and (v) the existence 
of complex spectral shapes, such as Figure 1 and Figs. 5 and 6 in Ryde et al. (2006). In 
addition, the thermal peak-energy is weakly dependent on the luminosity (T oc TL^^^) while 
the synchrotron peak-energy is much more sensitive to uncertain values of the magnetic field 
and of the energy in the electrons {Ep oc T^y'^^^B') and therefore expected to have a large 
dispersion. 

Indeed, from a theoretical point-of-view, thermal radiation originating from the base 
of the relativistic flow, where the densities are high enough to provide thermal equilibrium 
is an inevitable ingredient. This results from the fact that the fraction of explosion energy 
that is being converted to kinetic energy (that is later dissipated to produce the non-thermal 
spectrum) is necessarily less than 100%. Therefore, some fraction of the explosion energy is 
released in the form of radiation that thermalizes before escaping the flow at the photosphere. 
In addition, photons produced by any dissipation of the kinetic energy that occurs deep 
enough below the photosphere thermalize before escaping. The fraction of thermal emission 
that dissipative outflow gives rise to depends mainly on (i) the relation between the radius 
at which the flow saturates (thermal and kinetic energy are equal) and the photospheric 
radius, (ii) the efficiency of the non-thermal energy dissipation, and (iii) at what optical 
depth most of the dissipation occurs. In other words the existence of a thermal component, 
which is a natural outcome of relativistic outflow models, does not contradict the existence 
of a non-thermal emission, but adds to it. 

The interpretation of the prompt emission spectrum as being c ompos e d of a thermal 
component in addition to the non-thermal one, was put forward by iRyde I (120041 ). In this 
work, analysis of the time-resolved spectra of nine bright GRBs which were characterized by 



^We further discuss the difference and relation between our interpretation and the traditional Band model 
fits in S[43] 
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hard low-energy spectral slope, showed that a dominant thermal component could explain 
the observed spectra. The change in spectral shape over time, for instance by the low-energy 
spectral slope getting softer, was shown to be due to varying ratio of thermal to non-thermal 
flux and a change in the non-thermal spectral slope. This provides a natural explanation 
to the f act that the obse r ved s pectra are hardest during the initial phases of the pulse. 
Indeed, iGhirlanda et al. I (120031 ) already suggested that time-resolved spectra during the 
initial phases of some of the GRB pulses should be interpreted as pure blackbody emission. 

In iRyde 1 (2004) it was found that the observecEl temperature of the thermal component 
is approximately constant during the first few seconds, after which it decays as a power law 
in time T = T{t) oc t" ^, with powe r-law index ax ranging approximately from —0.6 to —1.1. 
It was later suggested (IRyde II2005I ) that the thermal emission component can be identified in 
many other bursts, in which it is not necessarily dominant over the non-thermal component. 



In this paper, we put forward the idea suggested by lRyde I (120051 ) . about the ubiquitous- 
ness of the thermal emission component in the prompt emission phase of GRBs. By using a 
large sample of 56 GRBs (see §2] for selection criteria) for which time resolved spectroscopy 
can be performed, we search for the existence of a thermal emission component and study its 
temporal behavior. In all the cases studied we were able to identify a thermal emission com- 
ponent. We extend here previous analyses of the blackbody (photospheric) emission. While 
former works were focused only on the temperature evolution, our present analysis includes 
the evolution of both the temperature and thermal energy flux, as well as the normaliza- 
tion of the blackbody fit, using a much larger sample than previously. This normalization 
corresponds to the ratio between the observed energy flux and the emergent flux from the 
photosphere, cxT'^, where a is Stefan-Boltzmann's constant. As we show below, it can be 
interpreted as the emitting surface, that is, the surface of the photosphere. We present the 
analysis method in §2 and summarize our results, in particular the spectral evolution of the 
blackbody component in §3. In §4, we discuss an interpretation of our results and elude on 
possible scenarios that could underly the observed properties. We summarize and conclude 
in 55. 



2. Analysis Method 



To properly analyze and determine the actual shape of GRB spectra one ideally needs 
burst detections with (i) high flux levels (ii) high time-resolution over individual pulses (iii) 



^In the following, all quantities are in the observer frame unless otherwise stated. E.g. quantities in the 
comoving frame are primed. 
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broad spectral coverage. Unfortunately, as of yet, all of these criteria are hard to meet at 
the same time, and the elucidation of the radiation process has not yet given a convincing 
answer. The best data for these investigations is still the BATSE sample from the Compton 
Gamma-Ray Observatory, which allows the analysis of spectra on short time scales over a 
spectral range of two order of magnitude above ~ 25 keV. Here, we focus on the study of 



3 ( 


Band et al. 


1993) 


Kaneko et al. 


2006) 



it is required to have a single or a couple of distinct pulses, which do not overlap each other 
significantly. Furthermore, we use a signal-to-noise ratio of more than 25 for the spectral 
analysis, which consequently determines the number and width of the individual time bins 
that are analyzed. The pulses are then required to have enough time bins for our temporal 
analysis to be meaningful, implying more than six time bins over the an alyzed duration. We 



also focus here only on the category of long bursts, longer than 2 seconds (IKouveliotou et al. 



19931 : iHorvath et al.ll2006l ). Our sample then consists of 49 pulses in 48 bursts (2 pulses are 
analyzed in BATSE trigger 2083). In addition to the complete sample of individual pulses, 
in order to demonstrate the robustness of our analysis method, we added to our sample eight 
bright bursts which have complex light curves with several heavily overlapping pulses. These 
are presented to illustrate the temporal behavior of the thermal component over the whole 
burst rather than only over a single pulse. The only requirement for these additional bursts 
is their brightness. The sample here is used for demonstration purposes, and a full sample 
of such bursts will be published elsewhere. In total we thus present the analysis of 56 long 
bursts, for which good time resolved spectra exist. The full sample is presented in Table 1. 

In the investigation below we follow the analysis by Ryde (2004, 2005) in which the 
background-subtracted photon spectra, N-E{E,t), are fitted by a Planck function combined 
with a power-law: 

where x = E/kT, kT is the color temperature of the blackbody, and k is Boltzmann's 
constant. The fitted parameters are thus the blackbody normalization A{t), temperature 
T(t), power-law index, s{t) and its normalization, B(t). The quality of the fit is given by a 
reduced x^- value. In this paper, we restrict our study to the thermal component, since it 
is the best constrained component and has a direct physical meaning (see §4 below). The 
non-thermal component is approximated by a single power-law over the analyzed energy 
range (25 -1900 keV). This assumption is, of course, too simplified. Nontheless it can be 
justified over the limited BATSE energy band. We further discuss this assumption in §4.51 
below; however a full, comprehensive study of the non-thermal component will be presented 
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in a future paper. 

The spectrum of the blackbody contribution to the observed flux is modeled with the first 
term in equation ([T]). At every instance (time bin) a temperature, T{t), and a normalization, 
A{t), are determined. The energy flux, F-BB{t) = J EN-^^dE, is given by an integration 
over all energies of this term after multiplying it by E{t). This integration yields -Fbb(^) = 
A{t) [kT^ vr^/lS. We next analyze the temporal evolution of the parameters T{t) and -Fbb(^)- 
We find that in 49 individual pulses it is possible to model the temporal evolution of the 
temperature with a broken power law, T oc t"^ before the break and T oc t^'^ after the break. 
We thus find it possible to extend the analysis by Ryde 2004 to a larger sample. As we 
show in §3, the temporal evolution of the flux can also be described by a broken power law, 
Fbb oc t"^ before the break and .F£3_oc t^f__after the break. Below we use the smoothly 



broken power law model derived in iRyde I (120041 ) and described by their equation (6) for the 



temporal evolution for both these parameters. The other bursts in our sample are complex 
and have heavily overlapping pulses and therefore an analysis of individual pulses is not 
possible. As described above for these bursts we make a similar analysis over the full burst 
duration. 

The data we analyze was taken by the Large Area Detectors (LADs) on BATSE. The 



data type we mainly use are the high-energy resolution burst (HERB) data (IFishman et al 



19891 ). These data have 128 energy channels over the energy range from 10 keV to a few 
MeV. The time resolution is multiples of 64 ms. The background estimates were made with 
the HER data, with a time resolution of 16-500 s. The spectral analysis was performed with 
the RMFIT package, version l.Obl (Malozzi et al. 2000) developed by the BATSE team in 
Huntsville, Alabama. The photon spectrum, Ne{E) was determined using a forward- folding 
technique; the spectral model was folded through the detector response matrix and was then 
fitted by minimizing the x^y using the non-linear Levenberg-Marquardt algorithm, between 
the model count spectrum and the observed count spectrum. This then gives the best-fit 
spectral parameters and the normalization. From these parameters the derived quantities, 
such as the energy spectrum and the energy flux were calculated. In the following the 
bursts will be identifie d by their BATSE trigger numbers according to the BATSE catalogue 



(jPaciesas et al. Ill999l ). see Table 1 
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Table 1. Temporal evolution of the thermal component in a sample of 56 bursts. Listed 
for the 49 studied pulses within these bursts are the power law indices (ap) and bx (bp) 
of the temporal evolution of the temperature (flux) before and after the break time to, as 
well as the power law index r of the temporal evolution of TZ. See the text for details. 



Burst 


Trigger 


a 7^ 




[s] 




bp 




r 


Comments 


910627 


451 


-0.2 ± 0.67 


-0.64 ± -0.32 


1.07 ± 1.3 


0.55 ± 0.36 


-2.67 ± 0.77 


1.73 ± 0.32 


0.78 ± 0.25 


a, pulse 2 


910807 


647 


-0.08 ± 0.05 


-0.77 ± 0.13 


2.95 ± 0.42 


0.5 ± 0.12 


-2.24 ± 0.33 


3.3 ± 0.2 


0.52 ± 0.05 




910814 


678 


0.001 ± 0.13 


-0.43 ± 0.043 


1.7 ± 0.6 


0.29 ± 0.28 


-1.43 ± 0.13 


1.67 ± 0.31 


0.02 ± 0.03 




910927 


829 














0.63 ± 0.05 


Complex; Fig. [9] 


911016 


907 


-0.18 ± 0.05 


-0.58 ± 0.02 


1.77 ± 0.19 


0.27 ± 0.12 


-1.26 ± 0.06 


1.82 ± 0.17 


0.54 ± 0.03 




911031 


973 


-0.19 ± 0.09 


-0.42 ± 0.06 


2.7 ± 1.2 


0.62 ± 0.12 


-1.66 ± 0.18 


2.85 ± 0.2 


0.54 ± 0.11 


a, pulse 1 


911118 


1085 














0.82 ± 0.03 


Complex; Fig. [9] 


920525 


1625 


0.4 ± 0.2 


-0.85 ± 0.1 


1.07 ± 0.12 


2.6 ± 0.6 


-3.5 ± 0.2 


1.27 ± 0.06 


0.39 ± 0.12 




920622 


1663 
















Complex; Fig. [9] 


920718 


1709 


± 0.5 


-1.2 ±1.1 


1.04 ± 0.9 


0.48 ± 0.9 


-2.2 ± 0.8 


0.69 ± 0.3 


0.38 ± 0.09 


pulse 2 


920830 


1883 




-0.43 ± 0.06 




± 0.1 


-2 ± 0.2 


1.45 ± 0.13 


0.71 ± 0.04 


b 


921003 


1974 


-0.05 ± 0.24 


-0.32 ± 0.03 


1.02 ± 0.62 


1.2 ± 0.2 


-0.8 ± 0.23 


1.33 ± 0.11 


0.61 ± 0.05 




921123 


2067 
















Complex; Fig. [TTI 


921207 


2083 


0.005 ± 0.07 


-1.01 ± 0.06 


0.977 ± 0.06 


1.23 ±0.1 


-2.8 ± 0.12 


1.12 ± 0.02 


0.76 ± 0.04 


pulse 1 


921207 


2083 


0.15 ± 0.09 


-0.77 ± 0.06 


1.68 ± 0.19 


0.98 ± 0.2 


-1.5 ± 0.2 


1.67 ± 0.17 


0.8 ± 0.11 


a, pulse 2 


930112 


2127 


0.03 ± 0.15 


-0.75 ± 0.16 


1.86 ± 0.42 


1.79 ± 0.43 


-3.1 ± 0.35 


1.77 ± 0.14 


0.71 ± 0.1 


a 


930120 


2138 


0.14 ± 0.05 


-0.43 ± 0.14 


1.98 ± 0.9 


-0.17 ± 0.8 


-1.8 ± 0.3 


2.3 ± 0.3 


0.14 ± 0.05 




930201 


2156 


-0.03 ± 0.12 


-0.4 ± 0.33 


4.24 ± 3.27 


1.75 ± 0.28 


-2.7 ± 0.2 


3.29 ± 0.18 


0.41 ± 0.07 




930214 


2193 


-0.25 ± 0.02 


-0.78 ± 0.04 


12.9 ±1.1 


0.15 ± 0.07 


-1.5 ± 0.1 


13.1 ± 0.8 


0.74 ± 0.04 




930612 


2387 


-0.09 ± 0.05 


-0.57 ± 0.06 


6 ± 1 


0.29 ± 0.07 


-2.1 ± 0.2 


8.8 ± 0.6 


0.43 ± 0.05 




940410 


2919 


-0.21 ± 0.2 


-0.69 ± 0.1 


0.85 ± 0.48 


-0.16 ± 0.24 


-1.13 ± 0.15 


0.89 ± 0.31 


0.44 ± 0.13 




940623 


3042 


-0.05 ± 0.06 


-1 ± 0.2 


5.3 ± 0.8 


0.78 ± 0.18 


-2.5 ± 0.3 


3.9 ± 0.2 


0.34 ± 0.11 




940708 


3067 


± 0.1 


-0.32 ± 0.06 


1.28 ± 0.46 


0.01 ± 0.11 


-2 ± 0.25 


3 ± 0.4 


0.07 ± 0.04 




941023 


3256 


-0.03 ± 0.16 


-0.66 ± 0.05 


1.54 ± 0.43 


0.05 ± 0.07 


-1.28 ± 0.08 


2.5 ± 0.2 


0.54 ± 0.09 




941026 


3257 


-0.18 ± 0.09 


-0.4 ± 0.03 


2.6 ± 1.1 


0.18 ± 0.1 


-1.34 ± 0.08 


4.38 ± 0.42 


0.48 ± 0.04 




950403 


3492 
















Complex, Figs. [Sl [Tol 


950624 


3648 


-0.04 ± 0.23 


-0.79 ± 0.07 


1.2 ± 0.3 


0.76 ± 0.93 


-1.83 ± 0.6 


1.41 ± 0.42 


0.7 ± 0.07 


pulse 3 


950701 


3658 


-0.27 ± 0.08 


-1.13 ± 0.15 


0.81 ± 0.11 


0.35 ± 0.2 


-4.5 ± 0.8 


0.95 ± 0.08 


0.72 ± 0.05 




950818 


3765 


0.13 ± 0.05 


-0.47 ± 0.05 


1.3 ± 0.15 


0.71 ± 0.09 


-2.35 ± 0.11 


1.48 ± 0.06 


0.02 ± 0.033 




951016 


3870 


0.02 ± 0.16 


-0.45 ± 0.07 


0.52 ± 0.19 


0.16 ± 0.25 


-1.5 ± 0.22 


0.77 ± 0.2 


0.22 ± 0.06 




951102 


3891 


0.52 ± 0.76 


-0.91 ± 0.5 


1.04 ± 0.47 


2.6 ± 1.6 


-2.4 ± 0.5 


0.95 ± 0.17 


0.41 ± 0.13 




951213 


3954 


0.1 ± 0.11 


-0.6 ± 0.1 


1.13 ± 0.33 


0.18 ± 0.18 


-1.37 ± 0.15 


1.1 ± 0.2 


0.01 ± 0.09 




951228 


4157 


0.07 ± 0.06 


—0.42 it 0.03 


2.5 ± 0.3 


0.76 ±0.1 


— 1.61 ± 0.11 


3.17 ± 0.14 


0.87 ± 0.08 


a, pulse 2 


960124 


4556 


0.07 ± 0.09 


-0.83 ± 0.15 


1.97 ± 0.3 


1.06 ± 0.2 


-2.59 ± 0.3 


2.01 ± 0.1 


0.39 ± 0.04 




960530 


5478 


-0.39 ± 0.07 


-0.68 ± 0.08 


3.4 ± 1.1 


0.09 ± 0.15 


-1.62 ± 0.15 


3 ± 0.4 


0.62 ± 0.04 




960605 


5486 


-0.1 ± 0.07 


-0.82 ± 0.1 


1.55 ± 0.26 


0.26 ± 0.18 


-1.63 ± 0.23 


1.66 ± 0.18 


0.56 ± 0.07 




960804 


5563 


0.04 ± 0.13 


-0.83 ± 0.21 


0.6 ± 0.12 


1.41 ± 0.25 


-3.8 ± 0.6 


0.7 ± 0.05 


0.69 ± 0.06 


a 


960912 


5601 


-0.17 ± 0.3 


-0.93 ± 0.09 


2.13 ± 0.6 


1.5 ± 1 


-1.7 ± 0.2 


1.61 ± 0.29 


0.76 ± 0.09 




960924 


5614 
















Complex, Figs. I8I10I 


961102 


5654 


-0.12 ± 0.06 


-0.62 ± 0.19 


18.1 ± 3.5 


0.78 ± 0.08 


-2.6 ± 0.31 


18.15 ± 0.79 


0.86 ± 0.08 


a 


970111 


5773 














0.38 ± 0.01 


Complex; Fig. [5] 


970223 


6100 


0.09 ± 0.36 


-0.46 ± 0.31 


1.18 ± 1.1 


0.62 ± 0.4 


-2.1 ± 0.3 


1.3 ± 0.26 


-0.02 ± 0.12 




970420 


6198 
















Complex; Fig. [5] 


970815 


6335 


-0.28 ± 0.02 


-0.56 ± 0.03 


3.51 ± 0.46 


-0.29 ± 0.09 


-1.38 ± 0.12 


3.04 ± 0.55 


0.39 ± 0.01 


subpulse at 3s 


970828 


6350 


0.19 ± 0.09 


-0.68 ± 0.13 


9.33 ± 1.14 


1.6 ± 0.3 


-1.93 ± 0.28 


9.32 ± 0.68 


0.44 ± 0.05 




970925 


6397 


-0.21 ± 0.04 


-0.55 ± 0.08 


6.4 ± 1.5 


0.37 ± 0.11 


-1.76 ± 0.11 


4.18 ± 0.32 


0.43 ± 0.07 


a 


971127 


6504 


-0.11 ± 0.09 


-0.67 ± 0.07 


3.8 ± 0.6 


0.21 ± 0.1 


-1.35 ± 0.08 


4.7 ± 0.4 


0.62 ± 0.03 




971208 


6526 
















very weak thermal flux 


980306 


6630 


-0.16 ± 0.04 


-1.14 ± 0.05 


1.58 ± 0.09 


0.31 ± 0.07 


-2.7 ± 0.2 


1.94 ± 0.07 


0.85 ± 0.06 




980718 


6930 


-0.3 ± 0.08 


-1.3 ± 0.2 


2.69 ± 0.2 


-0.12 ± 0.16 


-2.8 ± 0.2 


2.6 ± 0.2 


0.6 ± 0.1 




990102 


7293 


-0.36 ± 0.03 


-0.64 ± 0.04 


8.1 ± 1.3 


0.4 ± 0.14 


-1.78 ± 0.12 


5 ± 1 


0.52 ± 0.02 




990102 


7295 


-0.12 ± 0.06 


-0.74 ± 0.04 


3.3 ± 0.5 


0.83 ± 0.2 


-1.46 ± 0.15 


2.19 ± 0.2 


0.46 ± 0.07 


cxcl. subpulse at 10 s 


990123 


7343 


-0.09 ± 0.04 


-0.43 ± 0.15 


7.5 ± 2.2 


0.46 ± 0.33 


-1.15 ± 0.6 


5.94 ± 1.9 


0.25 ± 0.06 


pulse 1 


990316 


7475 




-0.315 ± 0.04 




0.71 ± 0.17 


-2 ± 0.22 


10.7 ± 0.9 


1.1 ± 0.22 


b 


990424 


7527 


-0.07 ± 0.03 


-0.76 ± 0.07 


2.5 ± 0.2 


0.08 ± 0.4 


-2 ± 0.6 


1.8 ± 0.7 


0.27 ± 0.06 




990816 


7711 


-0.41 ± 0.06 


-0.9 ± 0.07 


2.6 ± 0.5 


-0.01 ± 0.1 


-1.74 ± 0.14 


2.4 ± 0.2 


0.8 ± 0.03 




991216 


7906 


0.05 ± 0.05 


-0.8 ± 0.23 


0.98 ± 0.16 


0.94 ± 0.11 


-3.2 ± 0.7 


1.08 ± 0.08 


0.37 ± 0.04 


pulse 1 



"Break in TZ. 



No significant break in the fcT evolution 
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One can argue from a theoretical point of view that the spectrum we identify as a thermal 
component is formed either by photospheric emission without further dissipative effects (if 
scattering is the dominating opacity then the emitted emission is rather a Wien spectrum) , 
or by Comptonized thermal emission due to additional dissipation (IRees fc Meszaros II2005I ). 
Furthermore, the actual spectrum is most likely a convolution in time and space of several 
emitting regions, which, by necessity, broadens and dilutes the spectrum into a diluted 
and multicolor blackbody (see e.g. Pe'er 2008). Indeed, it is not statistically possible to 
dis tinguis h betw een a pure blackbody, a diluted blackbody, or a Wien spectrum as shown 
by iRyde I (120051 ). Therefore, it is still useful to fit this component by a Planck function, 
which in this sense captures such a peaked component, having in mind that this component 
probably is quasi-blackbody emission. 

Similarly, the single power-law used in the fits, approximates the accompanying broad- 
band, non-thermal emission. The limited band width prevents a detailed characterization 
of the non-thermal radiation, in particular its spectral shape. In a possible scenario, the 
photons emitted from the photosphere can serve as seed photons for C ompton scattering 



by ene r getic electrons pro duced by dissipation processes in the flow (e.g. iRees fc Meszaros 
( I2OO5I ): iPe'er et al. I (120051 )). If the dissipation occurs near or below the photosphere, energy 
exchange via both inverse and direct Compton scattering with the thermal photons may 
significantly modify the electrons energy distributi on, and as a cons equence a variety of 
complicated non-thermal spectra may be obtained (jPe'er et al. 1 120061 ). These can not be 



described by a simple power law over a broad energy range. However, a single power-law 
function is a relatively good approximation over a narrow energy band, such as the BATSE 
band. See further discussion the validity of this assumption in §4.5[ 



3. Spectral Evolution Analysis of the Prompt Phase Emission 

3.1. Decomposition into Spectral Components 

We show in Figure [1] the decomposition of the prompt emission of several BATSE bursts 
into the two spectral components, namely the thermal and the non-thermal components. 
These bursts have been chosen to illustrate the spectral dispersion that is observed. The 
thermal component is prom inent and is the cau se for the spectral peak, at an observed 



peak energy Ep ~ 250 keV (IKaneko et al. II2006I ). The broad band, non-thermal emission 



is captured by the power-law component. In this interpretation the spectral evolution is 
caused by the change in temperature, which can be significant both during the pulse as well 
as throughout the entire burst. We find that the temperature exhibits a characteristic decay 
behavior, a broken power law in time (see §3.2 below). This is similar to the behavior reported 
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by Ryde (2004). Different combinations of tlie tliermal and tlie non-tliermal components give 
rise to tlie variety of observed spectra; for instance, a variety of subpeak spectral slopes and 
peak energies. Spectral analyses, using the model in equation (1), on bursts observed with 
other satellite mission s (with rather limited band widths) have been reported in for instance. 



McBreen et all feOOGf ). Bosnjak et al. (2006), Falcone et al. (2007), and Bellm et al. (2008). 



Our fitting mo del is an alternatiy e to the Band model, which consists of two smoothly 



joined power-laws (IBand et al. 1119931 ). The two-component model gives fits with reduced 
X^- values around 1, which is similar to the xl values found for the Band model over the 
BATSE energy band (see Ryde 2004, 2005), however, having a somewhat larger dispersion 
in values. For the sample presented in Table 1, the xl values vary from 0.85 to 1.15 with 
typically 110 degrees of freedom per time bin. The averaged reduced for the whole sample 
is xl = 1.01. The lowest xl values (< 1) are obtained for time-resolved spectra in which 
there is emission at high energies that is not readily captured with a Band function. For 
these bursts the two-component model gives better fits to the spectra (see also Batellino et 
al. 2006 and Ryde et al. 2006). The highest xl values are caused by spectra which have large 
residuals at the lowest and the highest energies. In these cases the single power-law function 
does not fully capture the broad band emission at the extreme energies (see discussions in 
§4.51 for more details). 

Moreover, a majority of spectra fitted with the Band function in the BATSE energy 
band have f3 < —2. This is, for instance, the case for BATSE triggers 907 and 1663, shown 
in Figure 1. Our fits indicate, on the contrary, a rising high energy component. While 
these fits have similar as the Band function fits (in the BATSE range), their prediction 
of the flux above 2 MeV are quite different. One reason for the ambiguity of the exact 
spectral behavior in the BATSE range is the fact that at high energies the effective area of 
the detector decreases quickly. The rising spectral slope of course necessitates an additional 
high energy break. Broad band observations are needed to clarify this point. This will be 
further discussed in ^4.51 



The relative contribution of the thermal emission over the non-thermal emission is dif- 
ferent for different bursts. This is evident from the right-most column in Figure [H in which 
the temporal evolution of the ratio of the thermal to total emission is depicted. Here the 
thermal flux is the bolometric integrated blackbody flux and the non-thermal flux is inte- 
grated over the energy band ~25-1900 keV. The ratio found, ~ 30% — 50%, is thus an upper 
limit to the true (and unknown) bolometric ratio. It should be noted that while the relative 
strength of the thermal component varies in time, no strong recurring trend is revealed. 
For BATSE trigger 1663 the ratio remains approximately constant, while in other bursts 
it either increases or decreases. This is different from the temperature and flux behaviors 
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found below, for which a recurring behaviors are identified. On the other hand, these ob- 
servations indicate that, over an emission pulse, the thermal and non-thermal fluxes are not 
independent of each other, since the ratio remains roughly the same as compared to the flux 
variation themselves. 



3.2. Temporal Behavior of the Temperature 

The temperature of the quasi blackbody component shows a reccurring behavior. During 
the rise of the pulse, the temperature is typically constant or slowly declining, varying as a 
power law in time; the averaged values of the power law index is < ax >= —0.07 with 
a standard deviation 17(0^) = 0.19. After the pulse peak the temperature starts to decay 
faster, again as a power law; the average value of bx is < 6t >= —0.68 and ct(6t) = 0.24. 
This is consistent with the flndings in the works by Ryde (2004, 2005). Figure [2] shows eight 
bursts illustrating the variation that is observed. Some pulses have distinct breaks and the 
broken power-law model gives a univocal representation of the data, illustrated by triggers 
1709, 3042, 3765, 4556, 7527 in the flgure. Other bursts, however, have a smaller break, that 
is, the difference 6t — is small. An unusual behavior is exhibited by BATSE trigger 1883 
which does not have a measurable break in its decay at all. The measured values of and 
Bt, as well as the break time, to,T for the bursts in our sample are given in Table 1. 

The panels in the left-most column in Figure [3] show the distribution of the power 
law indices and bx for the sample studied here. The dispersions are rather large and 
for individual bursts these values can be rather different from each other. On an average, 
however, the temperature is constant until the break, at which it starts to decay as power 
law with an index, which is numerically close to —2/3 (see further discussions in §4 and also 
in Ryde, 2004). 

3.3. Temporal Variations of the Flux of Thermal Emission 

For each time interval during a burst, the temperature (T) and the blackbody normal- 
ization (A) are directly measured from a flt to the spectral data. The energy flux (Fbb) is 
subsequently calculated. Although the flux is not independent of the other two quantities, 
it has a physical meaning of its own and therefore, nonetheless the analysis of the temporal 
evolution of these three quantities (using minimization of the x^) ^'^'^ made independently 
of each other. 

During the temperature evolution described above, the energy flux contained in the 
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thermal component typically rises to the peak, which often coincides with the break in the 
temperature decay, after which it decays off. Figure H] shows the thermal fluxes for a number 
of pulses in our sample. As evident from the figure, the temporal evolution of thermal fluxes 
are described by power laws. The slopes vary from pulse to pulse, and during the rise phase 
the sample average of the power law index of < ap >= 0.63 with a standard deviation 
ai^ap) = 0.66 (following Fbb oc t^^^). During the pulse-decays, the corresponding average 
value of the indices are < bp >= —2.05 with standard deviation o"(6f) = 0.75. The broad 
dispersions around the averaged values are prominent in Figure [31 

As mentioned above, the break time in the temperature coincides with the peak in 
thermal flux (to,F) to within the errors. This is shown in the right-most panel in the lower 
row in Figure [31 The solid line shows a linear relation between the two times. 



3.4. Normalization of the Blackbody Spectrum 

In general, the parameter A{t) (see Eq. 1) describes the size of the photosphere from 
which the blackbody is emitted, as well as the efficiency of the emergent flux (see discussion 
below). Here, we deflne a dimensionless parameter TZ as ratio of the observed energy flux 
and the emergent flux (see also Pe'er et al. 2007) 




Therefore, the parameter TZ{t) is proportional to A(tY^'^. We now expand previous studies 
of the photospheric emission by studying the evolution of the parameter TZ(t) (instead of 
A{t)), since it contains important information regarding the emission site. 

The behavior of TZ over four individual pulses is shown in Figure [51 The insets show 
the count light curves over the same time intervals. A common property appears to be the 
power-law increase of TZ; TZ{t) oc Such a behavior is indeed the typical during most 
bursts, in particular, during their individual pulses. The ubiquitousness of this behavior is 
further shown in Figure [6l , [71 [HI For the bursts presented in Figure [6l the evolution of TZ are 
modeled by single power laws over the entire pulse duration. The value of (the dimensionless) 
TZ often varies over an order of magnitude, with typical values lying around 10~^^ (compare 
eq. [1]). As we will show in §4.2, this corresponds to a physical photospheric radius of 

fl^, = 2 X 10" ( ^ J (^) (^) f i±i) " cm (3) 

The characteristic temporal behavior of TZ shown in Figures [5l [6l [71 [S is very ubiquitous. 
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Since TZ is related to the photospheric radius, we find this behavior essential in understanding 
the properties of the photospheric emission, the emission site and the outfiow dynamics. 

In some bursts there is an indication of a break in the power-law increase of TZ after 
a few seconds, after which TZ is approximately time independent. Several examples are 
illustrated in Figure [3 Note that since the fiux levels become low at the end of the pulses, 
where these breaks in TZ occur, the error bars on the measured data points are large and 
consequently they exhibit large dispersion. Therefore, while the break is apparent, the exact 
behavior of TZ after the break is not very clear. Interestingly, we find that this break-time is 
independent of the break-time in temperature, and can be both leading (e.g. trigger 5654, 
with coincidence in trigger 2127) and trailing, with the latter being the typical behavior. 
The power-law indices prior to the break are on average comparable to those of the single 
power-law bursts shown above, indicating that the break might have a different origin. 

We have also identified a few bursts for which TZ is close to being time independent over 
the full duration of the pulse. The existence of such bursts is also evident from the histogram 
in Figure [31 where a few bursts show a power law index r close to 0. Several examples of 
these bursts are shown in Figure [HI For BATSE triggers 1625 and 5614, the constant level 
is reached after approximately 1 second (a few time bins); see further discussion in §4.2. 

In the investigation above we have studied individual puls es since these are the main con- 



stitue nts of a light curve and exhibit strong spectral evolution ( IBand et al. Ill993l : iRyde fc Svensson 



20021 ) . However, a large fraction of GRBs have spiky light curves in which individual spikes 
are too short to enable time resolved analysis. We have nontheless applied our method on 
a few bursts with highly variable light curves, bearing in mind that the time resolution of 
the individual time bins do not resolve the light curve variations. The inferred values of fiux 
and temperature in each time bin are thus averaged values, and as might be expected do 
not show any clear behaviors. However, we find that the temporal evolution of TZ{t) does 
show well defined behaviors. This is illustrated in Figure [9] where we show TZ = TZ{t) during 
six bursts of varying complexity, having several spikes which are not individually resolvable. 
The insets show the count light-curves over the same time intervals. A remarkable feature 
that can be seen in the figure, is that for many of these bursts, TZ{t) increases more or less 
independently of the light curve. This can, in particular, be seen for BATSE triggers 1085 
and 5773. Again a power-law increase in time is the prevailing behavior. Furthermore, for 
BATSE triggers 1663 and 6198, the parameter TZ{t) appears to asymptotically settle down 
to a constant value, during a period in which the light curve is in its most active phase. 
These findings demonstrate, first, that 7?. is a fundamental quantity and, second, that phys- 
ical information can indeed be extracted from bursts with highly variable light curves, by 
studying TZ. A more extensive study of complex hght curves will be presented elsewhere. 



-14- 



The histogram of the distribution of all the power-law indices r (from fits to 7^(t) oc t ^) is 
depicted in Figure [3] (upper, right-most panel). The distribution does not have a pronounced 
peak and the averaged value is at < r >= 0.51 while the standard deviation is a{r) = 0.25. 

The power-law fits to the 'JZ{t) data, both over individual pulses as well as complex 
bursts, often give good fits over the whole observed period (BATSE triggers 2193 and 7711 
in Fig. [5] and BATSE triggers 829 and 5773 in Fig. [9j). However, in some cases the power 
law only characterizes the overall (and indeed important) trend (e.g. BATSE trigger 1085 
in Fig. [H]). Variability on top of the power-law fit are apparent (correlated with the peaks 
in the light curve), but the magnitude of the variability is small relative to the range over 
which the fit is made. Similar deviations are seen for some of the bursts in Figure [6] (e.g. 
BATSE trigger 6630). More complicated fits (broken power laws etc.) could be made to 
analyze these secondary features. 

The investigation of the behavior of the parameter TZ above can be summarized as 
follows. Over an individual pulse structure, 71 typically increases as a power-law, extending 
over the entire pulse duration (exemplified by trigger 2193 in Fig. \5^. In a few bursts, at 
the end of the pulse duration, the parameter TZ settles to an approximately constant value 
(exemplified by trigger 3257 in Fig. [7]). This break time is independent of the break time in 
temperature evolution. For multi-pulse, and complex bursts we find, astonishingly enough, 
that the power-law can extend over the whole burst independent of individual flux variations 
(exemplified by BATSE trigger 5773 in Fig. [H]), while in bursts with well separated pulses 
the TZ power-law is clearly connected to the individual pulses (exemplified by trigger 2083 
in Fig. [6]). 

3.5. Comments on the Results 

The fits to the temperature, normalization (7^), and flux are made independently of 
each other. Since these parameters are related, consistency between the resulting power 
laws is expected. We indeed find that the power law indices for the flux, e.g. bp are close to 
the expected value, given the fits to T and TZ. However, by studying the power law indices 
of the decay phase of pulses in Table 1 it appears, at first glance, that for a few pulses 
the measured value of bp deviates somewhat from the expected value. The reason for this 
apparent inconsistency is often stochastical fluctuations in the flux data which leads to a 
measured break time which is somewhat different from that in the temperature evolution. 
This in turn leads to the deviation in the power law temporal slope. More importantly, 
though, is the late-time break sometimes observed in the TZ evolution. Consistency between 
the parameters should therefore only be expected if one restricts the fitted flux decay phase 
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to earlier times than the break time in 71. However, it is not a priori obvious from the flux 
curves that such a restriction should be made; the effect of the break in TZ on the flux curve 
is weak due to the large error bars on the data points. The conclusion is therefore that the fit 
results of the time evolution of the quantity Fbb, which is derived from the basic quantities 
A and T, can give spurious results in some individual bursts. 

As illustrated above, typically the flux peak is connected to the break in temperature 
(Fbb oc TZ'^T^) . In two cases, however, there is no observed break in the temperature (BATSE 
triggers 1883 and 7475). In these bursts, the flux peak is instead connected to the the break 
in TZ and is therefore different in origin. In Figure [7] the TZ evolution in trigger 1883 is also 
fitted with a broken power law with a break at the flux peak time. Such a fit is indeed 
consistent with the data. 

Finally, an inherent problem in describing data with a power law function is the sensi- 
tivity of the value of the index to the chosen zero point in time. In the fits made above the 
zero point was taken to be at the onset of the pulse. However, in many of the studied bursts 
the zero point is set by the trigger time of the detector, which does not necessarily exactly 
coincide with the onset of the pulse. This mainly affects the measured value of the pre-peak 
power laws, i.e. a-r and Oi?, which therefore should be treated with some caution. The post 
peak power laws are not affected by more than typically 5 %, by varying the chosen zero 
point around the onset of the pulse. 



4. Discussion 

As shown in this paper, we are able to identify the thermal emission component in 
all of the bursts studied so far, in the sense that we get acceptable fits to the data. The 
characteristic peak in the energy spectra of the prompt emission in GRBs is thus interpreted 
as in many cases stemming from a thermal component. Examples of such spectra are given 



in, for instance, iRyde I (120041 . l2005l ). In other cases, though, the power in the non-thermal 
components clearly has to peak at a higher energy, which is evident in BATSE triggers 1974 
and 7711 (Fig. [1]), and in trigger 7170 (Fig. [T2!) . For the latter cases, the thermal component 
contributes to the emission of the spectrum below E Fe peak. The most important result 
in our paper is that we find that the temporal behavior of the thermal component shows 
clear repetition in all the different bursts. Both the flux of the thermal component and 
its observed temperature show a broken power-law behavior in time: The temperature is 
observed to start off at a high level and thereafter it decreases, at first slowly, during the flux 
rise, and later breaking into a steeper decay, during the decay phase of the pulse. The break 
in the temperature occurs in conjunction with the flux peak. The parameter 7^, reflecting the 
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blackbody normalization, increases monotonically, typically as a power law in time. Most 
often, it has the same behavior over the entire pulse, and sometimes even over the whole 
burst. 



4.1. Consequences and Advantages of a Thermal Component Interpretation 



The existence of a photospheric component as 



deed expected both in dissipated kinetic outflows (IMeszaros. &: Rees 



part of a braod 



3and spectrum is in- 



2000a ) as well as in 



dissipative Poynting flux outflows (jPrenkhahnl 120021 : iGiannios fc SpruitI 120051 ) . The optical 
depth near the base of the flow is enormous, r > 10^^ (for a review, see, e.g., Piran 2005), 
and therefore photons emitted by the central engine or by any dissipation mechanism that 
occurs deep enough in the flow, necessarily thermalize before decoupling from the plasma 
at the photosphere. In addition, dissipation processes close to the photosphere (such as 
internal shocks or magnetic rec onnections) can enhance this thermal component by forming 
a Comptonized spectral peak (IRees fc Meszaros I l2005l ) . The original thermal peak can be 
upscattered up to a factor of 10 in energy due to the energetic electrons accelerated in the 
dissipation process. Further dissipation episodes expected during the prompt emission in 
the optically-thin region above the photosphere produce non-thermal photons. While the 
photosphere photons are, in principle, the first to reach the observer, in practice, due to 
Lorentz contraction the observed time-difference between these photons and non-thermal 
photons can be shorter than a millisecond and is thus not resolved. The observer is therefore 
expected to measure both components simultaneaously. The existence of a thermal spectral 
component in the prompt emission is thus inevitable in the standard fireball model scenario, 
and only its relative importance can be debated. 

The interpretation of the spectra as containing a thermal emission component in addi- 
tion to the non-thermal one, and the finding of the repetition in the temporal behavior of 
this thermal component have significant advantages over the purely non-thermal interpre- 
tation of the spectrum (the "Band function"). One obvious result is the ability to apply a 
consistent physical interpretation to the spectra in cases where the low-energy spectral slope 
is too steep to account for by the optically-thin synchrotron emission model. Another sig- 
nificant advantage of our interpretation of the data lies in the assumed physical origin of the 
thermal emission, as originating from the photosphere. For instance, the spectral peak has 
an immediate meaning as an effective (Lorentz-boosted) temperature. In the synchrotron 
interpretation of the Band model fits, the peak energy is instead related to the low-energy 
cut in the distribution of shock-accelerated electrons (characterized by an electron Lorentz 
factor, 7min) ; Ep oc T'j^^^B', where B' is the comoving magnetic field strength (neglecting 
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any dependency on the pitch angle distribution). There is therefore no a priori reason to 
expect clustering in Ep. While non-thermal emission originates from the dissipation of the 
flow kinetic energy, a mechanism whose details are unknown, much less uncertainty exist in 
the description of the photosphere. It is thus possible to study the properties of the photo- 
sphere, which is the innermost radius from which information can reach us, from study of 
the thermal emission. 

A correlation between the peak energy a n d flux naturally arises from the properties of 
the thermal emission (e.g. [Thompson et al. I (120071 )). Such correlations are indeed promi- 
nent within pulses and bursts. These are known as the hardness- int e nsity correlations ( or 
Golenetskii et al. (1983) relations; see also iBorgonovo fc Ryde I ( 120011 ): iButler et al.l (120071 )). 
The measured correlations in these works is though between the total flux (thermal -|- non- 
thermal) and the spectral peak energy, not only the thermal emission. However, in our 
interpretation these correlations originate from the underlying t hermal relatio n , wit h the 
additional non-thermal flux only affecting their appearance. The lAmati et al. I (120021 ) cor- 
relation represents a corresponding correlation between the flux and Ep values of the time- 
integrated spectra for an ensemble of bursts. The relation between temperature and thermal 
flux might be the underlying cause of this relation as well; see further discussion in §5 below. 

In many of the bursts studied here, the spectra are not dominated by the thermal 
emission, even within the observed energy band (25-1900 keV). We thus suggest, based on 
the results of this work, that thermal emission can in fact exist in a very large fraction, 
perhaps even all of the prompt emission spectra of long GRBs. 



4.2. Radiative Efficiency 



The radiative efficiency of the prompt phase emission is given by the parameter r] = 
Ej/{E^ + i^x) 5 where E^ is the radiated, prompt, gamma-ray energy (measured in the 
BATSE energy band) and i?K is the kinetic energy of the fireball right after the prompt 
phase. The estimation of Ey^ is somewhat dependent on the invoked afterglow model. How- 
ever, most estimates point toward very high values of the efficiency, ma i nly varying from a 
level of several percent to larger than 90 % (IFreedman fc Waxmarul200ll : IZhang et al.l 120071 : 
Nysewander et al.l 120081 ). This observational result seems to be in contradiction to the pre- 
diction of the internal shock model. In spite of its many successes, a major difficulty in 
this model is that, in general, internal shocks are inefficient in tapping the kinetic energy of 
the flow. Typically, only a few per ceiit of the fireball energy is converted into gam ma-rays 
JPaigne fc MochkovitchlEJ : iKumar 1^9991 : iBeloborodovlboOol : ISpada et al.lboooh . 
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Our new interpretation, and in particular the results presented here in Figure 1, may thus 
help to resolve this issue. According to this interpretation, thermal photons originate from 
the photosphere and thus they may originate, at least in part, directly from the explosion and 
not from dissipation of the kinetic energy in the flow. The fact that the thermal photons carry 
a significant fraction (30% — 50%) of the total observed flux in the BATSE band (see Fig. 1), 
therefore eases the energy requirements from kinetic dissipative models. This interpretation 
in supported b y num erical simulations of jet propagation out through the progenitor star by 



Lazzati et al.l (120091 ). They find that tangential coUimation shocks are important through 



out the jet, and that these generate continuous dissipation. This makes the jet internally 
hot producing a very bright photosphere. More than half of the total energy in the of the 
jet can thus be converted into radiation. 



4.3. Interpretation of the parameter 71 

The measured variations in the observed energy flux represents, to our interpretation, 
underlying variations in the energy input at the central engin^. The temperature, on the 
other hand, is given by the energy per particle and also informs us about the thermalization 
that takes place. While the flux and temperature have rather straight forward interpreta- 
tions, we will elude on the interpretation of TZ below. 

We assume that the flow is advected through a photosphere, at distance i?ph with bulk 
Lorentz factor F >> ©j^t ' where Gjet is the GRB jet opening angle. The observed flux, Fbb, 



is given by integrating the intensity over the emitting surface Fbb = (Stt/c^l) / djji iJiR^^J)^{aT' / 
Here T' is the comoving plasma temperature, which is related to the observed temperature 
via T = T' V^V = (F(l — /9/i))~^ is the Doppler factor, = arccos/^ is the angle to the line of 
sight, and P = {1 — F~^)~^/^ is the plasma expansion velocity. The ratio TZ = {Fbb / o'T^Y^'^ 
is calculated by integrating over 9: 

R^f^% (4) 

where z is the redshift and c^l is the luminosity distance. The c oefficient £ is a num erical 



factor of the order unity resulting from angluar integration (see iPe'er et al. I ( 120071 )). For 
bursts with known redshift z, the parameter TZ, which is a ratio of observed quantities, can 
thus be interpreted as an effective transverse size of the emitting region. 

According to this interpretation, a constant TZ means that the effective emitting area of 



■^as opposed to, e.g., obscuration. 
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the photosphere is time independent. A few of the bursts in our sample show this property 
and are thus the easiest to interpret. In Figure [8] several bursts for which this is the case are 
presented. We show in Figure [TD] the blackbody flux plotted versus its temperature for these 
bursts. These plots show the hardness-intensity correlations and, as expected, the observed 
flux is Fbb oc T^, with 6 ^ 4; the fitted values of the power-law indices are shown in the 
panels. This is the fundamental property of a blackbody emitter. The existence of this 
correlation illustrates directly the photospheric interpretation that is put forward here. 

We can further demonstrate how to extract information from bursts with complex light 
curves. Figure [TT] illustrates a specific example, namely GRB921123 (BATSE trigger 2067). 
The light curve of this burst consists of several overlapping pulses creating a complex light 
curve (see Fig. [H]). The tempo ral evolution of the temperature and flux track each other in 



this particular burst (see also ICrider et al. I (119971 )). The temporal evolution of parameter 



TZ is shown in Figure [91 From a statistically point-of-view the data can be fitted by a 
single power-law. However, another possible interpretation can be made. Instead of a single 
power-law, the TZ(t) data can be represented by two intervals of constant TZ values; one value 
during the interval between — 3 seconds and another value during the interval between 
6 — 15 seconds. Such an interpretation is indeed supported by the flux versus temperature 
plot (right-most panel in Fig. [TTl) . During both these periods the relation between the 
energy flux and the temperature is Fbb oc T*^ with g~ 4, which is indeed expected when TZ 
is constant. Thus in this interpretation, the effective transverse size of the emitting surface 
is constant during each time period, but changes between them. This burst illustrates how 
detailed analysis within the framework presented here can increase our understanding of the 
behavior of individual bursts. 



4.4. Temporal Evolution of the Thermal Component 

A major finding of this work is the well defined temporal evolution of both the tempera- 
ture and energy flux of the thermal component. The physical interpretation of this evolution 
is not obvious. It could be due to variation in the internal properties of the inner engine that 
produces the burst outflow (e.g. luminosity, baryon load, mass ejection rate, etc.), which are 
reflected in the observed temperature and thermal flux. 

A natural explanation to the late time temporal evolution (after the temperature break 
time) was suggested by Pe'er (2008). In this work, the properties of the photosphere in 
relativistically expanding plasma outflow, characterized by steady Lorentz factor F ^ 1 
were considered. It was shown there that the photospheric radius strongly depends on the 
angle to the line of sight, 6 (see also Abramowicz, Novikov & Paczyhski 1991). As a result. 
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thermal photons that decouple from the plasma at high angles to the line of sight, 6 ^ F^^, 
can be observed tens of seconds after the thermal photons originating on the line of sight 
(for typical parameters characterizing GRB outflows). Moreover, in this work the standard 
definition of the photosphere (as a surface in space from which the optical depth r = 1) 
was extended to consider the probability of thermal photons to undergo their last scattering 
event before reaching the observer from every point in space in which electrons exist. This 
enabled calculation (via probability density functions) of the late time thermal flux and 
temperature. The results found there, Fbb oc and T oc t^^/'^ — t^^/^ are remarkably 
close to the averaged values that we find for the fitted parameters (see Fig. [3]). The broad 
distribution of the TZ parameter found here may reflect the variation of the temperature 
evolution. 

The excellent agreement found between the theoretical prediction of the photospheric 
emission model of Pe'er (2008) and the data presented here, strengthen, to our view, our 
interpretation that indeed a thermal emission component exists during the prompt emission 
phase of GRBs, and that we were able to identify it correctly. The theoretical model, 
however, does not predict the early (before the break in the temperature temporal evolution) 
temporal evolution. The nearly flat behavior of the temperature at these times suggests that 
the emission is dominated by photons emitted on the line of sight, and that high- latitude 
emission effects can be neglected at early times. Thus, according to this interpretation, 
variation in the light curve during the early times directly reflect variations in the inner 
engine activity. This interpretation allowed us (Pe'er et al. 2007) to use data collected 
during this time to estimate the physical parameters of the GRB outflow, such as the size 
at the base of the flow and the outflow bulk Lorentz factor. 

We note though, that this interpretation does not naturally explain the lack of significant 
break in the temporal evolution of TZ at early times. Thus, alternative explanations may 
exist and further investigations are therefore ongoing. 



4.5. Broad Band Emission of the Non-Thermal Component 



While we have focused in this paper on the behavior of the thermal component, a full 
understanding of the energy release during the prompt phase in GRBs can only be achieved 
by a complete analysis of both the thermal and non-thermal components. In order to achieve 
this one needs a comprehensive analysis of the non-thermal component of time-resolved 
pulses over a broad energy range. While the observation s of the sub-MeV range is very well 



covered, on l y a few obseryation s exist at higher energies (lAtkins et al. Il2000l : iGonzalez et al. 



2003, 2009 



Abdo et al. Iboooh . Moreover, in most cases, these observations are integrated 
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over the entire duration of pulses and in some cases even over the entire duration of bursts. 
Since significant spectral evolution occurs during individual pulses, only indirect information 
regarding the emission processes can be found by studying such duration-integrated spectra. 
We therefore claim that the standard fitting of the GRB prompt emission, namely with the 
Band function, can only be regarded as a first approximation to the actual spectral form in 
the observed energy range. In addition as we argue here, for many bursts, the addition of 
a thermal component is required by the data, thus giving rise to two spectral breaks in the 
spectrum. 



4-5.1. Indications of multiple spectral breaks in the broad-band spectrum 



The BATSE spectra are often fitted with a Band function. One may thus argue that the 
broad-band spectral break in the prompt emission {E Fr) is at the break se en in the BATSE 
range, which on average is at approximately 250 keV (IKaneko et al. II2006I ). This is however 
not necessarily alway s the case. As shown in the spectral catalogue of BATSE bursts by 



Kaneko et al. I (120061 ) . approximately 10% of all bursts have, during their entire duration, 
time-resolved spectra with a high-energy power-law index j3 > —2, that is, increasing power 
with photon energy {E oc E^^"^). Moreover, approximately half of all bursts have periods 
throughout their durations over which (3 > —2. This means that the attributed peak energy, 
Ep, from the sub- Me V data is not the power peak of the broad-band spectrum, but rather a 
spectral break. A turn-over (peak) must exist at higher energies for these bursts. Similarly, 
the relatively low detection rate of burst emission in the Fermi LAT energy range, compared 
to what is expected from extrapolation of the BATSE /^-values, also indicate additional 
breaks in the super-MeV spectrum. The broad-band spectra are thus expected to typically 
have more spectral details than what can be captures by a Band function alone (see also e.g 
Barat et al. 1989). 



4.5.2. Broad-band spectral observations (keV-100 MeV) 



In a recent work, iGonzalez et al.l (120091 ) have completed the analysis of all 37 CGRO 
bursts for which both BATSE and TASC have a clear detection. Out of these bursts 9 have 
sufficient temporal resolution, giving a couple of spectra each (of a few second duration) 
during the rising phase of pulses. These are therefor e appropriate for direc t investigations of 
the underlying emission processes. In the work by I Gonzalez et al.l (120091 ) the spectra were 
fitted with a Band function and a smoothly broken power law. Here, we reanalyzed the 
BATSE data for these nine bursts over the same time intervals. We find that a thermal 
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component can indeed be identified in all cases. An example of our analysis is presented in 
Figure [H which shows the BATSE energy range (< 2 MeV) for GRB 981021 (trigger 7170). 
This presented spectrum is from an initial time bin, from 1 to 3 s after the trigger. We 
first present the Band function fit (Fig. [T2]left panel). In this fit t he high-energy spectra l 



index (3 = —2.87 and peak energy are frozen at the values found by iGonzalez et al.l (120091 ). 
Note that these values are obtained from their fit to both the BATSE and TASC data. The 
low-energy power-law index is a = —0.59 ± 0.01. The residuals between the model and the 
data, presented below the spectrum in the figure, do not show the random distribution that 
is expected from a model that captures the main physical properties of the data. Moreover, 
the p-value of the fit, with x^(dof) = 1.43(111), is p = 0.003. This clearly indicates the need 
for an extra component in the spectrum below its power peak at ~ 1600 keV. 

The right-most panel in Figure [T^ shows an alternative spectral fit to the same spectrum, 
using a thermal component and a power This fit results in kT = 160 ± 5 keV and a 

power-law index s = —1.17 ± 0.01, with = 1-03 (109). The distribution of the residuals 
has now improved, which is also reflected hj p = 0.38. Apart from improving the fit, the 
power law below the non-thermal peak at 1600 keV is now compatible with most non-thermal 
emission models, e.g. synchrotron emission. In the Band function fit the line-of-death for 
optically-thin synchrotron emission (slow cooling assumption) was violated {a > —2/3). In 
burst spectra similar to the one presented in Figure [T2| the non-thermal spectrum has a break 
at a few MeV, constituting the overall spectral power peak, while the thermal emission forms 
a shoulder on the low energy power law. Such spectra could be underlying the observations 
in e.g. GRB 900520a (Barat et al.l989) and GRB 080916c (Abdo et al. 2009), which exhibit 
spectral breaks at high energies. 



In some of the other bursts analyzed by iGonzalez et al.l (120091 ). the non-thermal com 



ponent does not have a detected break within the TASC energy range {E < 110 MeV). 
This is illustrated by the first time bin in GRB920311 (trigger 1473), for which the Band 
function fit to the BATSE and TASC data gives j3 = —1.71 ± 0.08. The spectral power thus 
increases with photon energy and exceeds that at the observed spectral break, E^^'^'^ = 390 
keV. Fitting the BATSE data alone, we find fcT = 68 ± 13 keV and s = -1.60 ± 0.05. The 
power-law index is consistent with the /3-value found over the TASC ranged. In this burst 
the spectral break is due to the thermal peak while the non-thermal component is consistent 
with a single power-law over the energy range 25 keV - 110 MeV. Similarly in GRB 930506 



a smoothly-broken power-law can be used instead, but since Ep is at 1600 keV, the difference is marginal 

^Indeed, forcing s — —1.71 for our fits over the BATSE range gives only marginally larger value and 
fcT = 83 ± 10 keV 
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(BATSE trigger 2329; Kaneko et al. 2008) there is no evidence for a power peak in the 
spectrum below 100 MeV. 

In many of the spectra taken from these nine bursts, our fits show that the thermal 
peak and the non-thermal spectral breaks are relatively close to each other. Nonetheless, we 
stress that these two peaks are not identical. For instance, in the interval between 1 and 3 
seconds of GRB 920311 (1663), we find kT = 73 keV and = 715 keV. The model 

used to fit this spectrum consists of a Planck function and a smoothly-broken power-law. 
An indication of such a high-energy break is indeed evident in the fits of this burst in Figure 
1. In these bursts the spectral power in the non-thermal component is lower than in the 
thermal component. As a result it is easy to misinterpret the spectrum as consisting of only 
a Band function, with a peak at the thermal peak (see Fig. [T]). 

As mentioned above, the spectra from the nine bursts studied here are from the initial, 
rising phase of the pulses. As shown by Ryde (2004) the non-thermal power-law typically 
softens during the duration of a pulse (see also Fig. [1]) , and therefore a different broad-band 
spectral behavior is expected at later times of individual pulses. However, there are no 
available data for these nine bursts, which sufficiently resolve the decay phase of the pulses. 

We point out that the two spectral breaks identified and discussed here, one due to 
kT and the other one due to E™"^'^^', will by necessity be obscured if the interval studied 
are time-integrated over significant spectral evolution. This is because both peaks vary in 
time. Moreover, obliging eff ects of the derived photon and energy spectra further obscures 
the apparent spectral shape (IFenimore et al.lll983l : iBromm fc Schaeferlll999l ). The fitting 



is namely performed in count space and the photon flux data points, used in the plots, are 
estimated based on the best fit model. Different models therefore lead to slightly different 
derived photon fluxes (note however that the spectral fitting is unaffected by this obliging 
effect). This underlines the need for time- resolved and forward folding analysis for a proper 
understanding of the emission processes during the prompt phase in GRBs. 

Even though only limited amount of data are available for time-resolved, broad-band 
spectra, and these are mainly for hard bursts (having large Ep values), we draw an important 
conclusion. While the dispersion in the energy of the thermal peak is relatively small, the 
peak energy of the accompanying non-thermal component has a much larger dispersion, 
ranging from ~ 700 keV to larger than 110 MeV. 
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4-5.3. Clues to the origin of the nonthermal component 



We thus conclude that our interpretation, which is an alternative to the Band function, 
is consistent with the observed data. In this work we do not attempt to fully analyze the 
non-thermal part of the spectrum (this will be presented elsewhere). However, we briefly 
review some suggestions on the origin of the non-thermal component. In the dissipative 
outflow model, the non-thermal component is expected to have a non-trivial, broad-band 
shape: The energy dissipation above the photosphere can lead to accelerated particles 
(emitting synchrotron and/or invers e Compton emission) t hrough Fermi (stochastic) ac- 
celeration, or streami ng instabilities (IFrederiksen et al.l 120041 ) and/or converter acceleration 
( iDerishev et al.ll2008l ) in the flow. Plasma oscillations induced through Compton scatteriri g 
by the burst emission itself, can also give rise to accelerated electrons (iFrederiksenI |2008| ) . 
These lead to a variety of spectral slopes and several possible spectral breaks. 

Furthermore, depending on the optical depth that the photospheric photons experience 
at the dissipation site the particle distribution will be reprocessed by the photons, thereby al- 
tering the emitted spectrum to a varying degree, partly thermalizing the spectrum: Thermal 
photons could serve as seed photons for IC scattering by energetic electrons which are pro- 
duced by the dissipation mechanism that dissipates the outflow kinetic energy. In this case, 
a variety of non-thermal spectra may be produced, of which the exact shape depends on the 



prop e rties of the unspecified dissipation mechanism (IRees fc Meszaros I l2005l : iPe'er et al. 



20051 . 12OO6I ). In particular, the thermal component may or may n ot dominate the to tal ob- 
served flux, depending on the values of the free model parameters ( IPe'er et al. II2006I ). Such 
reprocessing of the non-thermal particle distribution by the photosphere photons introduces a 
relation between the two components. In addition, since the non-thermal emission is drawn 
from the kinetic energy of the outflow it should be correlated with the thermal emission, 
apart from efficiency factors which depend on how and where the dissipations occurs. 



Alternatively, if the outflow is Poynting-flux dominated (lDrenkhahnll2002l : lGiannios fc Spruit 



20051 ) t he dissipation occur s at rn agne tic reconnection sites throughout the flow. In such a 



scenario iGiannios fc SpruitI (120071 ) and lGianniod (120081 ) have shown how the fraction of ther- 
mal emission in the spectrum varies with the Lorentz factor of the outflow. For moderate 
Lorentz factors a thermal component can claim a substantial fraction (~ 30%) while for 
very large Lorentz factors the fraction is low. To summarize, in order to understand the 
largely unknown dissipation processes giving rise to the non-thermal component, temporally 
resolved emission pulses, observed over a broad energy range, are required. 
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4.6. Further Analyses 

Here we have studied the spectral evolution of a complete sample of long pulses in GRB 
light curves. We have also studied a sample of 8 bursts with complex light curves and showed 
in §4.31 and in Figure [9] that in spite of that the pulses overlap each other, we were able to 
extract information about the photosphere. Therefore a more exhaustive study of bursts 
with complex light curves is expected to reveal further information. 

On a few occasions, the Swift satellite has provided a broad spectrum (0.2- 150 keV) of 



the prompt phase, with its two instruments BAT and XRT (IGehrels et al. II2004J ). However, 
to extract physical information on the photosphere these detections must be sufficiently 
strong to enable time-resolved spectroscopy in a satisfactory manner. Moreover, a wider 
spectral range is required. The Fermi Space Telescope covers an energy range of approxi- 
mately 10 - 200 GeV with its two instruments (Burst Monitor and the Large Area Telescope) 
and thus is expected to provide a sufficiently broad energy range (see Band et al., in prep.) 
to allow this analysis. 

In addition, the late-time fla res, which are obser ved in the afterglow, can be further 



studied with this model (see, e.g. iFalcone et al.l (120071 )). Once available, using the method 



described here and in Pe'er et al. (2007), the data will allow to discriminate between com- 
peting flaring models, such as slower propagating shells or late-time engine activity. 

The strong photospheric emission identified here can be accompanied by a unique neu- 
trino signal. This takes place if two conditions are met: first, the dissipation of the fireball 
kinetic energy occurs below or near the photosphere (see Pe'er, Meszaros & Rees, 2005, 
2006), and second, the dissipation process produces population of energetic protons. The 
energetic protons then interact with the low energy thermal photons via photo-meson and 
Bethe-Heitler interactions, as well as proton-proton interactions (see Koers, Pe'er, & Wijers 
2006) to produce high energy neutrino signal. The fiux of energetic neutrinos that is pro- 
duced by these processes depends on the spectral shape of the interacting photons. Thus, 
a thermal photon spectrum gives rise to a unique neutrino signal (see Wang & Dai 2008, 
Murase 2008). The results obtained in these works are different than former estimates (Wax- 
man & Bahcall 1997, Dermer & Atoyan 2003), which considered pure power law spectrum 
of the interacting photons, neglecting thermal component. A full analysis of the neutrino 
signal expected during the prompt emission phase has however to consider the full spectral 
shape of the interacting electrons, hence requires first a full modeling of the prompt emission 
spectrum (both thermal and non-thermal components). 
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5. Conclusions 

We have identified a thermal component in the prompt spectra of GRBs, with a temper- 
ature of a few 100 keV. This component has a characteristic behavior, consistent with that 
of the flow photosphere. This emission is superimposed on a broad-band, non-thermal emis- 
sion. The broad band spectrum therefore consists of at least two peaks, or spectral breaks, 
which are more or less pronounced, one from the photosphere and one from the non-thermal 
emission. The power peak of the of the broad band X-ray and 7-ray spectrum can either be 
that of the photosphere or that of the non-thermal component, for strong bursts lying at a 
few MeV. 

We have shown on a sample of 56 long BATSE bursts (25-1900 keV), that it is possi- 
ble to model instantaneous spectra of the prompt emission with a two-component model, 
consisting of a Planck function combined with a power law, approximating the non-thermal 
component in the analyzed energy range. Our interpretation implies that the peak in the 
photon spectrum, at ~ 250 keV, is attributed to the quasi blackbody emission from the 
photosphere of the relativistic outflow. Theoretical arguments that lie in the heart of the 
GRB fireball model easily imply that a blackbody component is expected. Furthermore, the 
values of the physical parameters derived from our analysis (e.g. T, Fr r, -Rph, F; see also 



Pe'er et al. 2007) are consistent with the predictions of the fireball model (IRees fc Meszaros 



2OOOI ). In summary, the arguments for this interpretation are the following: 



1. Goodness of fits. From a statistical point-of-view the two-component model fits the 
data well which is reflected by the reduced xl — ^ values. In some cases our model gives 
better fits than the commonly used broken power law (Band et al. 1993) model (see §3). In 
other cases the Band model gives a better x^-value, indicating the need for a more complex 
model (other than a single power law) to fit the non-thermal component. 

2. Radiative Efficiency. While the non-thermal photons originate from dissipation of 
the kinetic energy, e.g., by internal shock waves, or magnetic reconnection, thermal photons 
can originate directly from the progenitor, i.e., photons that are released in the explosion 
and thermalize before decoupling at the photosphere. As we have shown in Figure 1, thermal 
photons carry 30% — 50% of the total flux. Thus, inclusion of the thermal photons in the 
total energy budget of the prompt emission can contribute to explaining the high efficiency 
reported, which is difficult to account for by theoretical models of kinetic energy dissipation. 

3. Recurring behavior. The behavior of the thermal emission component is similar for 
most bursts, showing a particular and recurring behavior. For the 49 pulses we studied 
we found that both the temperature, T(t), and the flux of the thermal component, -Fbb(^), 
exhibit a well defined, and most importantly repetitive behaviour: a broken power law 



-27- 



in time. During the pre-peak phase the temperature is approximately constant while the 
energy flux rises. During the post peak-phase, both the flux and the temperature decay 
as power-laws in time; Fbb t''^ , and T oc t''^. While the power-law indices do have 
a broad distribution as shown in Figure [3] and Table 1, the sample- averaged value of the 
temperature indices is < br >= —0.68 and the average of the flux indices is < bp >= —2.05. 
The repetitive behaviour found, together with the excellent agreement between the results 
found here and the theoretical predictions of Pe'er (2008), strongly support, to our view, the 
interpretation of the prompt emission spectra as containing a thermal component. 

4. Blackbody normalization. The normalization of the thermal emission (parameterized 
by 7^(t) = {Fbb/o'T^Y^'^) shows a distinct, also recurring, temporal behavior which is largely 
independent of the temperature evolution. We found that the parameter TZit) typically 
increases in size as a power law in time over the whole pulse and sometimes even over the 
entire duration of bursts with complex light curves. We have thoroughly discussed ( §3.41 §4.3[ 
Figs. El [6], O [8], [9]) the evolution of TZit) and have shown that this parameter, which is directly 
derived from observed quantities, is directly related to the photospheric radius. Hence, it 
carries with it a significant information on the physics of the prompt emission. In principle, 
the photosphere is the innermost radius from which any photon can reach the observer; 
therefore, studying emission from the photosphere holds the key to our understanding of the 
prompt emission and the properties of the site from which it emanates. 

5. Spectral shapes. Our interpretation of the prompt emission alleviates many of the 
problems of purely non-thermal interpretations, such as the optically-thin synchrotron inter- 
pretation. In particular, hard sub-peak spectra that challenge the optically-thin synchrotron 
interpretation of the prompt emission spectra are naturally obtained in our interpretation. 

6. Peak energy. Identification of the sub-MeV spectral peak energy with the thermal 
component ~ kT, implies that it is independent of the local magnetic field strength and par- 
ticle acceleration mechanisms, in contrast to the non-thermal (synchrotron) interpretation. 
This facilitates explanations of observed peak energy clustering and spectral evolution. In 
the analysis carried out in §??, we show that the break in the non-thermal component, ac- 
companying the thermal emission, is observed to have a larger dispersion as may be expected 
by theory. 

7. Spectral correlations. Thermal emission gives a natural correlation between tempera- 
ture and flux. This correlation should be underlying the hardness-intensity correlation (HIC 
or Golenetskii relation) observed within bursts, as well as the Amati-like relations which are 
found for ensembles of bursts and therefore relate them to each other. We note that for 
several of the studied bursts in this work the observed flux varies as the fourth power of 
the observed temperature, which is expected from a blackbody emitter. The implication for 
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the Amati-like relations, however, should be taken with some care: First, these relations are 
derived for time-integrated spectra, which is different from the time-resolved spectra con- 
sidered here. And second, we expect Componization of the thermal component to alterate 
the position of the peak (Pe'er et al., 2006). Quantitative studies of this alteration require 
knowledge of the distribution of the non-thermal electrons produced by the dissipation pro- 
cess close to the photosphere. This, in turn, can be modeled by study of the non-thermal 
part of the spectra. A detailed study is on-going. 

We conclude that the characteristic behaviors that exist for the temperature and the TZ 
evolutions arc the defining properties of the photosphere, and must therefore hold the key 
to our understanding of the prompt emission and the properties of the site from which it 
emanates. 
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Fig. 1. — Examples of time resolved E Fe spectra, observed by BATSE, fitted with the two- 
component model; solid line is the Planck function and the dashed line is the non-thermal 
(power law) model. The small panels beneath each spectrum show the residuals in units 
of the standard deviation, a. The spectra are rebinned to have a signal-to-noise ratio of 
unity to make the plots clearer. The full resolution is however kept for the residuals panels. 
For every burst, which is identified by its BATSE trigger number, two time-bins are shown 
in the first two columns. BATSE trigger 1663 has a complex light curve, while the other 
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Fig. 2. — Examples of temperature evolution from individual pulses in different bursts. A 
distinct break is here seen in all cases but trigger 1883 which is best fitted by a single power- 
law. The slopes vary somewhat from pulse to pulse, but the general behavior is universal. 
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Fig. 3. — Parameter distributions of the analyzed bursts. The panels in the left-most column 
depicts the histogram of the power-law indices, before (a^) and after (67-) the break in the 
temperature curve. The column in the middle correspondingly depicts the indices before 
(ap) and after (fep) the flux peak. The upper right-most panel depicts the histogram of the 
power law index of the TZ evolution. Finally, the lower right-most panel shows the correlation 
between the temperature break time and the flux peak time. The red line shows a linear 
relation. See the text for further details. 
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Fig. 4. — Examples of blackbody flux evolution from individual pulses in different bursts. 
A distinct peak is seen in all cases. Both the rise and decay phases are describes by power 
laws. Note that the actual rise index is very dependent on the zero point of the abscissa. 
In many cases the trigger is such that data is only available after part of the flux rise has 
occurred. The red dots refer to additional pulses that are excluded in the presented fits. 
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Fig. 5. — Evolution of the parameter TZ = (Fbb/'^T^Y^'^ describing the ratio between ob- 
served, thermal flux and the emergent flux. The corresponding count light curves are shown 
as insets (arbitrary units). A remarkable power law is exhibited, much independent of the 
rise and decay of the pulse. 
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Fig. 6. — More examples of the evolution of the parameter TZ = {Fbb / ciT^Y^"^ describing 
the ratio between observed flux and the emergent flux for individual pulses. A universal 
behavior is seen in all bursts. They exhibit a power-law increase over the whole pulse. 
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Fig. 7. — Examples of the evolution of the parameter TZ for burst with an apparent change 
in the power law index r. These bursts initially exhibit a power-law increase, similar to the 
bursts shown in Figs. [5] and [H In these bursts, after a few seconds the increase ceases and 
TZ becomes close to constant. 
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Fig. 8. — Examples of the evolution of the parameter TZ in bursts for which r is close to zero. 
These bursts are indeed consistent of having TZ being constant throughout the analyzed time 
period. The data points marked by red dots in triggers 1625 and 5614 are probably part of a 
short initial i ncrease. The light curves of triggers 3492 and 5614 reveal a complex structure 
(see Fig. 6 in lRyde fc SvenssonI (120021 )). 



-42- 



lE-18 - 



lE-19 




lE-18 



lE-19 





lE-18 



6E-19 



4E-19 



2E-19 



10 



8E-19 



6E-19 



4E-19 



\ 5773 


if J 






-r 1 

■ 




A 


1 


10 20 
, . , , 1 




2E-19 - 



9 10 11 12 



Time [s] 



Fig. 9. — Evolution of the parameter 7^ = [Fbb / (^T^y/'^ for bursts with complex and heavily 
overlapping light curves. The corresponding count light curves are shown as insets (arbitrary 
units). A remarkable power law is exhibited in many of the bursts, much independent of the 
complexity of the light curve. 
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Fig. 10. — Plots of the observed blackbody flux, Fbb, as a function of the temperature, 
kT, for the bursts pre sented in Figure El Such c orrelations are also known as the hardness- 
intensity correlations (jBorgonovo fc Ryde II2OOII ). The energy flux is found to increase to the 
fourth power of the temperature, and is therefore proportional to the emergent flux, crT^. 
The values of the individual power law fits are shown in the panels. The existence of such 
bursts illustrate directly the photospheric interpretation of the spectral peak. 
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Fig. 11. — Spectral evolution of the blackbody emission in BATSE trigger 2067. Its complex 
light curve is shown in Figure M (a) Thermal flux light curve, (b) Temperature evolution, 
(c) Flux versus temperature. Further discussion is given in the text. 
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Fig. 12. — -E-Fe spectrum from an initial time bin (1-3 s) in GRB981021 (BATSE trigger 
7071). The panel to the left shows the spectrum fitted to the Band fun ction with [3 set to -2.87 
and E-p = 1591 keV found from the TASC fit (iGonzalez et al.ll2009l ). The residuals clearly 
indicate the need for an extra component in the spectrum {x^ = 1.43 (111)). Moreover, 
the fitted value of a = —0.59 ± 0.01 violates the synchrotron line-of-death (a > —2/3). 
The panel to the right shows the spectral fit with kT = 160 ± 5 keV and power law index 
s = — 1.17±0.01, giving = 1-03 (109). The fit is improved and the non-thermal component 
now has a < —2/3, which is consistent with synchrotron emission. The change in flux scale is 
due to the obliging effect of the deduced photon spectra (see the text for further discussion). 



